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Abstract 

The comprehensive mass insertion formula of the b — > sqq transitions was 
provided for the chargino contributions. The method clarifies physics behind the 
calculations of the super symmetric contributions to FCNCs. The result was applied 
to the analysis of the difference of the CP asymmetries of the B —> Kit decays. It 
is shown that the supersymmetry is not responsible for the discrepancy. 



1 Introduction 



The CKM picture of the standard model (SM) has been confirmed by experimental and 
theoretical developments [H |2j [3j H] . Although a lot of the experimental results of flavor- 
changing neutral currents (FCNCs) are consistent with the SM, there are processes whose 
measurements deviate from the SM prediction [5]. Some of them have been reported in 
the b — s transition processes. Definitely, it is too early to conclude the anomalies due 
to physics beyond the SM unless uncertainties especially from hadronic contributions are 
reduced. Nonetheless, it is important to explore the new physics (NP) in the flavor sector, 
because the FCNCs are sensitive to flavor violations in TeV scale. 

Supersymmetry (SUSY) is one of the most attractive extensions of the SM. It involves 
new flavor structures due to scalar partners of the SM matters. They contribute to the b—s 
transitions by gluino and chargino diagrams [6 J . There are two methods to describe them: 
i) the loop diagrams are evaluated without approximations [7] , and ii) they are expanded 
by small parameters (see, e.g., [B]). The former provides full evaluations, while the latter 
has virtues that the contributions can be understood clearly. In fact, physics behind 
the former calculations becomes clearer, and cancellations which appear in the former 
method such as the GIM mechanism are performed automatically. In the literature, the 
latter description, called the mass insertion (MI) approximation, has been discussed for 
the gluino contributions. The MI formula for the chargino contributions to b — > sll is 
also given in Ref. J8J, though some typos were recently reported [9J. That for b — > sqq are 
found in Ref. [TUJ ITT] . However, some contributions are not included, as will be shown 
later. In this paper, the comprehensive chargino contributions to the Wilson coefficients 
of b — > sqq will be provided in the MI approximation. 

The analysis will be applied to the Kir puzzle. The B — > Ktc decay is dominated by 
the penguin contribution and sensitive to the NP (see, e.g., [12] for a review). Currently, 
the measurement of the CP asymmetry of the B° — > K + n~ decay is different from that of 
B — > K^tt at more than the 5a level pQ. In the SM, such a large deviation is unexpected 
(see [12J for a summary of the theoretical predictions), which is called the Kn puzzle. 
In the topological decomposition of the amplitudes [13], the difference originates in the 
color-suppressed tree (C) and electroweak penguin (P ew ) contributions. The experimental 
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results may indicate enhancement of C and/or P ew by NP contributions. 

The SUSY solution to the puzzle by means of the chargino contributions has been 
discussed in [TJ1 [151 IB] • In contrast to the gluino contributions which have already been 
constrained [TTj, the authors insisted that the chargino can induce a large CP asymmetry 
of B — > K7r. However, it was reanalyzed in [18] with the full loop calculations, i.e., i) in 
the above category, and found to be less significant. This conclusion will be confirmed by 
the MI approach. Also, the contributions with the squark mixing of (5^)32 are analyzed, 
which were discussed in [HI (T5J [16]. They can be included by turning on the charm 
Yukawa coupling and will be shown to be too small to solve the K7r puzzle explicitly by 
using the MI formula. 

This paper is organized as follows. The Wilson coefficients of the chargino contribu- 
tions as well as those of the gluino are provided in Sec. [2j They are applied to the CP 
asymmetries of B — > Kir in Sec. EJ The results are summarized in Sec. HI 

2 Wilson Coeffients 

The b — > sqq transitions are represented by the following effective Hamiltonian [19], 
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^eff - -7=- 2^ \ 



V2 

p=u,c 



CiO? + J2 c *°* + c Tr°Ty + c *gOsg + (L R) 



+ h.c, (1) 



_i=l,2 i=3 

by ta „ g 0Ut heav, d e,ees of fr ee dom , where f „ th e . 0r d qUar , He r e, K . 



defined by the CKM matrix as X p = V P bV*. The operators are defined as| 



q=u,d q=u,d 
q=u,d q=u,d 

q=u,d q=u,d 



1 There arc additional operators in the effective Hamiltonian up to the dimension six such as those 
with the scalar or tensor Lorentz structure. In this paper, we are interested in the processes that are not 
suppressed in the chiral limit. Since the extra operators involve the chirality flip of the quarks, they are 
irrelevant, and the operators in ([5]) are enough. 



3 



g=u,d 



q=u,d 



o 



n 



16tt 2 



o 



9, 



8G 



16tt 2 



(2) 



where the color indices are denoted by a and 0, and the electromagnetic and strong 
coupling constants are e and g s , respectively. In the electroweak (EW) penguins, C i=7 - W , 
e q is the electromagnetic charge of the quark, q. The above operators are often called the 
left-handed currents, while the chirality-flipped operators and their Wilson coefficients, 
Oi and Ci, are called the right-handed currents. 

Focusing on the transitions between the second and third generations, there are eight 
types of the squark mixings, which are classified by the chirality index of the squarks 
and whether the mixing belongs to the up- or down-type squarks. Let us consider the 
superCKM basis, where the quark supermultiplets are rotated such that the quark mass 
matrices are diagonalized. Off-diagonal components of the squark mass matrices, — £ = 
b*(Am 2 ~)ijSj + t*(Aml)ijCj, are parametrized as (5^)32 = (Am?) i: ,/m?, where i,j denote 
the chirality of the squark, L and R, and m~ is a diagonal component of the squark mass 
matrix. Thus, the SUSY contributions are represented by the eight MI parameters, 



32- 



(032, (S d LR ) 32- 



(Si 



LL/32, 



RL) 32, 
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In this paper, the soft squark masses are simply set to be same except for the right-handed 
stop, and rn| is chosen to be this common mass. For the down-type squark mixings, it is 
convenient to use their Hermitian conjugation, (6^)23 = (^J^- 

When the down- type squark mixings are present, the gluino-squark diagrams dominate 
the SUSY contributions to b — > sqq. The Wilson coefficients are obtained as [61 [20] : 
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where Bx^Px^, and M 0) f, are the loop functions (see App. [A} with x = m|/m~. Here 
and hereafter, the Yukawa coupling of the strange quark is neglected. The contributions to 
the right-handed currents are obtained by flipping the chirality such as (8l L ) 2 3 — > ($101)23 
and ( 5^)23 — > (8rl)23- ^ is noticed that the magnetic dipole terms, CV 7 and Cgc have 
two types of the contributions for (Sf jL )23 are (5 RR ) 2 3. It is easy to see that the tan/3- 
enhanced one is usually dominant, which is often called the double MI contribution or 
the induced LR mixing. 

For the four-quark operators, the gluon box and the gluon penguin contributions are 
included in (J3J). There are additionally the EW penguins by mediating the photon 
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while those to Cs and C10 are small. The contributions to (7, are obtained by flipping the 
chirality. The following Z-penguin contributions are also added to the coefficients, 
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The loop functions, Z a and Z , are found in App. [A] The Z penguins with (5 ( [ L )23 and 
(fiim) 23 can De comparable to (JS]) when tan/3 is large, though they are proportional to 
mj. Lastly, the box diagram contributions to C7 and C$ can be finite for (5j jL )23 and 
(5^)23, when there is a large difference between the sup and sdown masses. Since they 



are smaller than 



by 



[m 



m| )/m? they are discarded in this paper. 



The up-type squark mixings appear in the chargino contributions to b — > sqq. The 
charginos consist of the Wino and Higgsino, which mix to each other by the electroweak 
symmetry breaking. Since the weak boson masses are much smaller than the SUSY scale, 
the chargino amplitudes can be expanded with respect to mw,z /w so ft as well as (fifj)s2- 
The Wilson coefficients of the four-quark operators are classified as 
1 
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where C 9 represents the gluon penguins, C 7 the photon penguins, Cz the Z penguins, and 
Cb u and Cb a the box diagrams. In the MI approximation, they become 



c, 



c> 



c 



L>u 



c, 



M 2 W 



27ix\ t m\ 



1 L 



(fi U LL)^ L (x W ) 



+ (fiR L )s2 



fm t M 2 „ RL1 

g \ X W> -V> x t R ) 1 ^2 



m tA t p RL2 



( x Wi x t R ) 



Or 



M 2 W 



2lTx\ t ml 



q l 



(fi U LL )32P y LL (x W ) 



+ (fiR 



RL) 32 



a,. 



V 



X 



m- 



w ■ 



x in %t R ) H "7 ' 



4rr A t sin 



(*2 



m t M 2 m t A 



m 2 - 



m- 



pLL t _ 
r Z \ X W> A" tt 



+ (fi U R 



RL) 32 



07 



sin m= 



(fit 



LL)32 



m t A 
m\ 



Pz ( x w> x t R ) 



Mj LL1( , m t M 2 m t A t LL2 



m- 



m , 



m- 



+ (fiR 



RL) 32 



ft, 



\mtM 2 RLl 

\ x Wi "V' x t R ) * 



M 2 ™t A t dRL2 



2 2 



-^u ( x Wi X t R ) 



sin" 771= 



(fit 



LL)32 



B d \ x w) H s — I 



6 



RL) 32 



\™ t M 2 RL1 . 



1 



t B RL2 (r~ t- 1 



(7) 



where the trilinear coupling of the scalar top, A t , is defined as £ = —m t A t i R t L + h.c.. The 
contributions to Cg and Cio are suppressed. The chargino contributions to the magnetic 
dipole operators are obtained as 
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All the loop functions are summarized in App. [Aj and the parameters are defined as 
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where M 2 and /i# are the Wino and the Higgsino mass parameters, respectively. In the 
chargino contributions, the right-handed stop mass, m tfl , is left explicit, because the Z- 
penguin contribution becomes larger when it is hierarchically small. Since the strange 
Yukawa coupling is neglected, the right-handed currents do not receive corrections. 

In the superCKM basis, the CKM matrix appears in the chargino-quark-squark 
vertices, which is another source of the flavor violation. However, this effect is discarded 
in this paper, because we are interested in the cases when the corrections are dominated 
by the squark mixings. 

There are no chargino contributions from {5 RR )3 2 and (5^)32, when the Yukawa 
coupling constants of the second generation are neglected. This is because the Wino 
couples only to the left-handed quarks. When the charm Yukawa coupling is turned on, 
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the contributions become finite. The leading contributions are obtained as 
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where the loop functions are in App. [A] The other contributions are negligible. Although 
they are suppressed by the charm Yukawa coupling, since the constraints become weaker, 
the mixings can be sizable. It is found that the contributions to C*7 7 and Csg are not 
suppressed by M^/m?, but rather, they include tan/3. 

The Wilson coefficients are set at the scale where the heavy fields decouple, i.e., the 
weak scale for the SM and the SUSY scale for the SUSY contributions. They evolve from 
the input scale down to the m& one by solving the renormalization group equations (see, 
e.g., [19J). Then, the b — s transition amplitudes are evaluated. In the next section, the 
gluino and chargino contributions are studied on the CP asymmetries of the B — > Kir 
decays in the MI approximation. 

Before proceeding, let us compare the MI formula to the full evaluation of the loop 
diagrams. The full method means that the loop diagrams are evaluated in the mass 
eigenstate basis. The calculations are easy to preform, whereas it is relatively hard to 
understand the underlying physics. This is because the full formula includes cancellations, 
e.g., by the GIM mechanism. In contrast, the MI formula is suited for understanding the 
process, because the mass matrices are expanded, and such cancellations are absent. Since 
the calculations are involved because there are a lot of diagrams after the mass matrices 
are expanded, we compared the results of the MI formula to those of the full method 
in order to confirm them. The full evaluations are found in the literature (see [18] as 
a recent one). We checked that there is a good agreement between the two approaches 
in most of the parameter space. Exceptionally, when the right-handed stop is much 
lighter than the other superparticles, a subleading contribution can be dominant, because 
a corresponding loop function is enhanced. Then, the next-to-leading order evaluation is 
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required. Otherwise, the MI approximation works well. 

In the MI approximation, a part of the chargino contributions has been missed in the 
literature [10]. In particular, the contribution with A t has not been found, though those 
to C*7 7 and C$g are sizable. In the SUSY models, A t is considered to be comparable to 
the gluino mass, for instance, by the renormalization group evolution. Further, the recent 
LHC searches for the Higgs boson [21] may indicate relatively large A t [22]. Therefore, the 
contributions could be crucial for studies on the FCNCs. Moreover, the chargino mass 
matrix has not been expanded in [10J. Since the mixing between the gaugino and the 
Higgsino stems from the electroweak symmetry breaking, the off-diagonal components of 
the mass matrix can be expanded, which is considered in the present result. This enables 
us to avoid the cancellations which appear in [TO] . 

3 B — » Ktv 

The difference of the CP asymmetries of the B — > Kir decays has been measured as pQ 

AA CP = A CP (B + K + n°) - A CP (B° K + tt~) = 0.124 ± 0.022, (10) 

where the errors are added in quadrature. This is unexpected in the SM. The amplitude 
of the each decay mode can be decomposed by the topology of the weak transitions 



-A(B° K-n+) ~ P + 2 -P^ + T (11) 



P 



2 

1 + _ r g we ^Ew + rre ^T e -^3 



-V2A(B- -> K-n°) ~ P + P EW + Hp£ w + T + C 

2 

1 + r EW e iSEW + -rg w e i5 Ew + (r T e il5T + r c e i&c ) e"^ 3 
3 

where the QCD penguin P, the color-allowed EW penguin Pew, and the color-suppressed 
EW penguin P^ w are approximately proportional to A c in the SM, while the color-allowed 
tree T and the color-suppressed tree C are proportional to X u . Hence, a relative weak 
phase, X u /\ c oc e - ^ 3 , arises at the second line in the each amplitude. The strong phase 
relative to P is denoted by 5j. The ratios of the magnitudes are defined by 

r E w = \Pew/P\, rg w = |PeV/P|, r T = \P T /P\, r c = \P C /P\. (12) 
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In the decay amplitudes, the hadronic matrix elements are evaluated by using the QCD 
factorization as [231 [18] , 

P ~ X c (a c 4 + r^a c 6 )A nK , T ~ \ u aiA %K , C ~ \ u a 2 A Knj 

3 3 
Pew ^ (~ ai + ° 9 ) ^ K7r ' ^ew - ( a io + r f a §) (13) 

where subleading contributions including the weak annihilation are omitted for simplicity. 
The definitions of the parameters are found in [23]. The SM naively satisfies 1 > Tt, tew > 
r EW' r c ( see > e -g-> [II])- The difference of the CP asymmetries is represented as 

A^4 CP ~ — 2r c sin 5 C sin (fi 3 + 2r EW r T sin(5 EW + 8 T ) sin 3 . (14) 

Since the SM predicts r T , r EW = O(0. 1) and r^ w , r c = 0(0.01), Av4 CP ~ 0.01 is expected, 
which is smaller than the experimental result, fllOp . 

If the NP is responsible for the discrepancy, the SM evaluation ffl^j) is affected. The 
NP contributions are parameterized as 

P+lpg w + T^P + ^pC w + T + X m , 

-Pew + 0—7- Pew + C + Y^p, (15) 

where Xnp represents the NP contribution common to B° — > K~tt + and B~ — > K~tt°, 
while Inp contributes only to P~ — > K~tt°. Since the NP generally involves new CP- 
violating phases, there are four types of the corrections to AAqp'- 

1. When Xnp is dominant with a sizable weak phase relative to P(SM), the isospin 
violation is from P E w = Pew(SM), and AAcp is evaluated as 

AAcp(NP) ~ 2r EW r x sin(<5 EW + <5x)sin0 x , (16) 

where r x ,8x and (fix are |X NP /P(SM)|, the strong and the weak phases relative to 
that of P(SM), respectively. Thus, Xnp ~ P(SM) is required to account for ffTQ]) . 

2. When X^p is large, but its weak phase is almost aligned to P(SM), Xnp contributes 
to P as P — > P + Xnp. Then, AAcp has the same form as f lT4|) . In order to enhance 
it, Tj in ( 1T4|) needs to be magnified. Thus, X NP ~ — P(SM) is required. Note that 
the branching ratios of B — > Ktc are also suppressed. 
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3. If Ynp dominates the NP contributions with a large weak phase, AAqp becomes 

Av4 C p(NP) ~ -2r y sin 5 Y sin <p Y , (17) 

where the parameters are defined similarly to ffl6l) . Then, |1np| is required to be at 
least comparable to |Pew(SM)| for ffTUj) . 

4. If Fnp is large, but the weak phase is suppressed, Ynp contributes to Pew- Noting 
that the CP violation originates in T(SM), AAqp becomes 

AA C p(NP) ~ 2r Y r T sm(5 Y + <5 r )sin0 3 . (18) 

The experimental result ({TO]) indicates that Ynp is as large as P(SM). 

It is commented that when the NP appears in the right-handed currents, their contri- 
butions are evaluated by extending Ci in [23] as C; —> Ci — Ci, where the relative sign 
is determined by the parity of the final state [21]. Although the SUSY is unlikely to 
contribute to the strong phases [25], since there are discussions [18] . we focus on the 
magnitude of the SUSY contributions in the following study. If none of the conditions, 
[THU is satisfied, the measured difference ( |T0l) is not explained by SUSY, irrespective of 
the strong phase. 

In many cases, the SUSY contributions are dominated by the chromomagnetic dipole 
term, C^g- According to the MI formula, it is larger than the four-quark operators. From 
the relation between a« and Ci |23j, P in (fl3|) receives a correction as 

«4 + rfa 6 = -^{Ik + r«)Cf G + ..., (19) 

where Ik is J cfe $x/(l — x), and $^ is the leading-twist distribution amplitude of the 
K meson. It is emphasized that this contribution corresponds to X NP , whose weak phase 
is from the squark mixings in Csg- 

The SUSY contributions to C 8 g is tightly constrained by Br {b — > sj). This is because 
the structure of C*7 7 is very similar to Csg, as shown in Sec. |2] explicitly. Since the 
experimental results [1] agree well with the SM prediction [25], extra corrections are 
restricted to be within the range, 

-0.28 x 1(T 4 < AB(6 37) < 1.08 x 10" 4 , (20) 
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600 800 1000 1200 1400 

m soft [GeV] 

Figure 1: Contours of |Xnp/P(SM)| are shown for 0.2, 0.4, 0.6, 0.8 and 1. The gray 
shaded region is excluded by Br(6 — > s-y) at the 2a level. The parameter m so f t represents 
the soft breaking parameters and fin- The squark mixing other than Im {5f iR )23 are set 
to be zero. 



at the 2a level. This means that CV 7 is around 0.3 at the scale. In addition, the recent 
measurements of B — > K*fi + n~ support that CV 7 Cg is consistent with the SM, where Cg 
is the Wilson coefficient of the effective operator, Og = (sl7/^l)(^7^) [27]. Since the 
sign of Cg is not flipped by the SUSY contributions [28], CV 7 is dominated by the smB 

In practice, when (<5^ L ) 23 and (5^)23 are considered, CV 7 and C 8 g are larger than the 
four-quark operators by tan/3 (see (j4])). They contribute to X NP and b — > 57. In Fig. [TJ 
|^np/-P(SM)| and the b — > 37 bound ( 120]) are shown for m so f t and ^(5^)23 x tan/3. It 
is found that Xnp is less than 30% of P(SM) under the constraint, and the maximum is 
insensitive to m so f t . On the other hand, Ynp is found to be tiny. If the photon penguins 
are compared to C^g, © is proportional to a e without tan/3. Although the Z penguins 
can be comparable to the photon penguin, they do not dominate the EW penguins at least 



2 The amplitude also receives contributions from the chargino and charged Higgs diagrams with the 
flavor violation by the CKM matrix. They are assumed to be cancelled to each other for simplicity. Even 
if they are included, the following conclusion does not change qualitatively. 
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for tan/3 < 50, because they are proportional to m h (see (jSJ)). Hence, Y N p is estimated as 

Y NP ~ P E w(SUSY) ~ - 27TNc ^ P(SUSY) ~ -^-P(SUSY), (21) 

CfOs ( I K + r£ ) c; 8G tan (5 

where the loop functions are taken into account in the last relation. Since tew(SM) is 
of order 0.1, Y NP is found to be less than percent of Pew(SM). Therefore, none of the 
conditions, HHH is satisfied, and the measurement ( ITUj) is not explained, because of the 
tight constraint by b — > sj. 

The situation is very similar for (5^)23 and (5^)23- The magnetic dipole terms in 
(J4]) are enhanced by rrig/mb, whereas Inp from ([6]) is proportional to tan f3 without 
l/(Gi?m|). Then, Cj^/Csg in (J2TJ) is estimated to be ~ G^m^ tan (3, which leads to 
Ynp/Pew(SM) < 1% under the constraint from b — > sj. This is too small to explain f fTUj) . 

In the case of the chargino contributions with (<52l)32> the magnetic dipole terms 
in ([8]) include tan/3, though they are proportional to M^. As seen from (J7J), the Z 
penguin is not suppressed by M^, whereas it is proportional to a e w? t . Thus, C 7 $/C%g ~ 
a e mf / (M^ tan (3) is obtained for (l2Tj) . The SUSY contributions are controlled by CV 7 
and Csg, and f lTOj) is not explained by {S^ L )^2 by b — )■ 37. 

The b — > 57 constraint is relaxed, when (8^)32 is considered and if A t is suppressed, 
as found from ([7]) and (]8]). There is no enhancement for C$g compared to the gluon 
penguin, C g , in the SUSY contributions to the decay amplitudes. Since the Z penguin 
is not proportional to M^, the ratio in f l2~Tj) becomes C 7)9 /C 8 g ~ a e m?/M^. Thus, Y NP 
can be larger than Y NP ~ P(SUSY) for rriq ~ ITeV, though it decreases as rriq is lowered. 
In Fig. |2j contours of |Ynp/Pew(SM)| and 6—7-57 are plotted for m so ft and Im(<5^ L )32. 
It is seen that Ynp can reach about 40% of Pew(SM), namely a few percents of P(SM), 
for m soft ~ 150GeV. Then, AA CP (SUSY) could be as large as 0(0.01) from (JT7D, if the 
b — > S7 bound is imposed. 

There are constraints on ( 5^)32 other than from 6—7-57 such as B — > K*fi + fi~ [9] 
and the vacuum stability [29]. According to the recent analysis based on the updated 
experimental data of B — > K*fi + fi~ [9J, (5^)32 is limited to be less than ~ 0.1 for 
^-soft ~ lOOGeV and ~ 0.3 for m soft ~ ITeV. The vacuum stability condition [29J provides 
(3rl)32 ^ 0.5 and 0.3 for m soft ~ 500GeV and ITeV, respectively. Thus, Ynp is bounded 
to be less than 10% of Pew(SM) according to Fig. [2j Thus, (5^)32 cannot be large enough 
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100 200 300 400 500 

m so ft [GeV] 

Figure 2: Contours of |>Np/P ew (SM)| are shown for 0.1, 0.2, 0.3 and 0.4. It does not reach 
0.5. The gray shaded region is excluded by Br(6 — > sj) at the 2a level. The parameter 
m soft represents the soft breaking parameters and while A t = is taken. The squark 
mixing other than Im {5r L )32 are set to be zero. 

to explain ( flOl) . Furthermore, since direct searches for the superparticles already excluded 
the light mass region [30] . the bound is likely to be much severer, and it is concluded that 
(^rl)s2 cannot account for f fTOj) . 

In the above study on (5^)32, A t was assumed to vanish. As the coupling increases, 
the magnetic dipole operators are enhanced (see OH])). Then, the b — > S7 bound becomes 
severer, and Fnp is disallowed to be large. On the other hand, when the right-handed 
stop is much lighter than the other superparticles, the loop function of the Z penguin 
increases (see App. |A|. However, the bound from B — > K*{i + fx~ simultaneously becomes 
severer [9]. We could not find the region where the difference of the CP asymmetries of 
the B — > Ktc decays ( TTUj) is explained by (5^)32- 

Let us comment on the right-handed mixings of the up- type squarks, (S^ R ) 3 2 and 
(5^)32. As found in (Q, their contributions are proportional to the charm mass, m c . 
Since the contribution to CV 7 and C$g are suppressed, the mixings can be as large as 
unity. However, it is estimated that Y NP reaches less than 1% of Pew(SM) even for 
r^soft = lOOGeV, tan/3 = 40 and Im^^,^ = 1. Thus, it is impossible to solve the Kn 
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puzzle by (5% R ) 32 and 

In the QCD factorization, one of the largest uncertainties stems from parameterizations 
of the divergences in the soft gluon interaction and the weak annihilation (23] . So far, the 
scenario S4 in [23] was applied to the estimation of ( TTSj) . If the parameters are varied, 
A^4cp can change by ~ 5% (see, e.g., [H]). Nonetheless, we checked that the SUSY 
contributions are irrelevant for the puzzle. 

The above study is compared to the previous works which investigated the chargino 
contributions for the Kit puzzle [TU [151 HSJ HE]- The authors in [HJ [T5j [16] insisted 
that the puzzle was solved, whereas those in [18] obtained the opposite conclusion. The 
difference exists in the method to evaluate the loop diagrams: the former used the MI 
approximation, while the latter relied on the full method. In this paper, we revisited the 
Kit puzzle with the MI formula and obtained the same conclusion as [18]. Although the 
MI formula used in [TJJ [15] missed some contributions which mentioned in the previous 
section, they may not be the leading reason of the contradiction. We found that some 
numerical results in [Tl} [15] are not reproduced, which could cause the difference. On the 
other hand, the strong phase could be a source of the gap between [TU [TSJ [T6] and [IB] , as 
mentioned in [18]. In the former analysis, the phase is taken to be free, while it is fixed in 
the QCD factorization in the latter. In this paper, the puzzle was studied with the strong 
phase of the SUSY contributions supposed to be a free parameter. It was concluded that 
the chargino contributions are tightly restricted irrespective to the strong phase. Thus, 
it is believed that the SUSY is not responsible for the puzzle. 

Let us comment on a loophole of the above study. The SUSY contributions could be 
cancelled to each other for Cj T Then, large squark mixings might be considered to explain 
the anomaly. In that case, careful analyses are required, because such large mixings easily 
spoil other experimental constraints. For instance, when (S'[ L )23 is as large as unity with 
(Sf jR )23 tuned to cancel CV 7 , SUSY contributions to B s — > J/4>4> as well as Am s may 
conflict with the experimental results [5]. Anyway, the complete study is devoted for 
future works. 
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4 Conclusion 



The MI formula of b — > sqq was provided for the chargino contributions as well as those of 
the gluino. This approach has an advantage that they are easier to understand compared 
to the full evaluation. For instance, the cancellations which appear in the latter method 
such as the GIM mechanism are performed automatically, and it is possible to find out 
relevant contributions. 

The MI formula was, then, applied to the analysis of the Kir puzzle. It was found that 
the SUSY models are not appropriate to explain the experimental result of the discrepancy 
between the CP asymmetries of B° — > K + tt~ and B — > K ± tt°. Since the structure of 
the SUSY contributions is clearly seen by means of the MI approximation, the conclusion 
is expected to be less dependent on details of the models. 

Let us comment on a future prospect of the Kir puzzle. The CP asymmetries of the 
B — > Kir decays follows a sum rule unless the anomaly is due to the EW penguin |31j . 
If the violation of the sum rule will be observed in future, physics beyond the SM is 
indicated other than the SUSY models. 
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A Loop Functions 

In this appendix, the loop functions in the gluino and chargino contributions are provided. 
A.l Gluino contributions 

The loop functions that appear in (TJ]), ([5]) and are given by 

1 + Ax — 5x 2 + Ax In x + 2x 2 In x 



B 1 (x) 
B 2 (x) 



8(1 - xf 
5 — Ax — x 2 + 2 In x + 4x In x 
2(1 - a;) 4 ' 
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Pi(x) 
P*{x) 

M 2 (x 
M 3 {x 

M 4 (x 
M a {x 
M b (x 

Z a (x) 

Z b (x) 



1 - 6x + 18x 2 - 10x 3 - 3x 4 + 12x 3 In x 
18(x- l) 5 ' 

7 - 18a: + 9x 2 + 2x 3 + 3 In x - 9x 2 In x 
9(x- l) 5 ' 

= 4^ (a), 

-1 + 9a; + 9a; 2 - 17a; 3 + 18a; 2 In x + 6a; 3 In x 
: I2(x- l) 5 ' 

— 1 — 9a; + 9x 2 + a; 3 — 6a; — 6a; 2 In a; 
1 6(x- l) 5 ' 

= -3M 4 (x), 

3-3x 2 + (1 + 4x + x 2 )lnx 





(x - l) 5 


2x 3 + 3a; 2 


— 6x + 1 — 6x 2 In x 




6(1 -x) 4 


x 4 — 8x 3 - 


1- 8x - 1 + 12x 2 logx 


12(1 -x) 5 



A. 2 Chargino contributions 
A.2.1 The functions fW, g\ j \ hi 

The chargino contributions are represented by the functions, f- j \ g!f\ hf 



n 



(0 



(0 



(0 



(0 



(0 



AO 



n 



(i 



(i 



ff\x,y) = -2yf™(x,y), = iTM), 



[x,y) 



1 



x(logx — logy) 



x-y (x-y) 2 

(x,y) = fi°\y,x), / 2 (0) (x) = / 2 (0) (x,l), 

( r v ) - x + y x(\ogx-\ogy) (0) _ (0) 

{X ' V) ~ 2y(x-y) 2 (x - yf ' h [X) ~ h ^' ij ' 

(x,y) = fi°\y,x), /f(x) = / 4 (0) (x,l), 
1 + x 



;x) = (^iy + (^T^ logx ' 

5 + x 1 + 2x , 

+ 7 — logx, 



(x) = 

(x,y) 



2(x-l) 3 (x-1) 4 
3x — y x 2 (log x — log y) 



, A (1) (x) = / 1 (1) (x,l), 



2(x — y) 2 (x — y) 3 

(x,y) = fi%,x), / 2 {1) (x) = / 2 (1) (x,l), 
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(1 



(1 



If 
ff 



ft 



(2 



(2 



(2 



The definit 



[x,y) 



2x 2 + 5xy — y 2 x 2 (log x — logy) 



6y(x - yf (x - y) 4 

(x,y) = fi 1 \y,x), / 5 (1) (x) = / 5 (1) (x,l), 
. . 1 + 5x x(2 + x). 

(x) = -2(^T)i + (^Tji l0 ^' 

(*) = -if 

l + 10x + x 2 2(1 + a) 



+ 



logx, 



[x,y) = 



3x(x - l) 5 (x - l) 5 
llx 2 - 7xy + 2y 2 x 3 (logx - logy) 



6(x — y) 3 (x — y) 4 

(x,y) = / 1 (2) (y,x), / 2 {2) (x) = / 2 (2) (x,l), 



(2), 



(x) = 



1 - 5x + 13x 2 + 3x 3 x 3 
12(x- l) 4 (x - l) 1 

-17 - 8x + x 2 



■ log X, 



6(x- l) 4 



1 + 3x 
(x - l) 5 



on of the functions g 



9 S 0) 


(^,y,^) 






(^,y,^) 






(^,y) = 


9 r 


si" 


(x,y) = 


si" 


9?' 


(x,y) = 





x logx 



(x-y)(x-z) 
x 2 log X 



(x - l)(x - y)(x - z) 



logx. 

+ (x ^ y) + (x ^ z), 

+ (x <->■ y) + (x <H- z), 



The definition of the functions h,,: 



hi(x,y) 


= fi°\x,y)-f[ 1 \x,y), 


hi(x) 


= hi(x,l), 


h 2 (x,y) 


= fi 1 \x,y)-fi 2 \x,y), 


h 2 (x) 


= h 2 (x,l), 


h(x,y) 


= / 2 (1) (x,y)-/ 2 (2) (x,y), 


h 3 (x) 


= h 3 (x,l), 


h\{x) = 


2/ 3 (0) (x)-3/| 1) (x), 






h 5 (x) = 


3/ 4 (1) (x)-4/ 3 (2) (x), 






h 6 (x) = 


fi 1) (x)-ff\x), 






h 7 (x) = 


xf£\x)- 2 -ff\x), 
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2 

h 8 (x,y) = xff\x,y) - -/ 2 {2) (x, y), h 8 (x) = h 8 (x,l). 



A. 2. 2 Penguin 

The loop functions for C g , C 7 and Cz in (J7|): 
P?t*) = -\f?l*\ 

p rli (t „ , _ 9 f?\x)-f?\x,z)-(x^y) 
9 M '^" 8 (a._ y )(i_ z ) 



9 K ,yj 8{l-y)\ 3 V ; l-y 
P 1 LL (x) = -Af! i 2) (x)- 9 -h 7 (x), 

-'-7 \X,y,z) 



(x-y)(l-z) 
9 f /2 (1) (») - /2 (1) ( x , z) g[ 1} (x,z)- (x,y,z) 



2(x — y) \ 1 — z x — y 



^r(x,y) = - Utf>(x 



l-y I ° v ' i-y 



2(1-1/) I 1-Z/ 

pLL( „ x _ -fl^fo g) - (x^y) 

(x-y)(l -*) 
, /s (&) - 9i (x, z) g { 2 ] (x, z) - g)p (y, z) 



2(x-y)(l-z) 2(x-y) 



2 



pmif \ _ 9? {x, z) - g { 2 ] {y, z) + y{ 1} {x,z)~ y[ 1} (x, y, z) 
yx,y,z)— — r 

2(x - y) 

,ffl2 L „a 1 ( Al)f„\ „(!)/„. „A , „f(0)/N l,(i). 



The loop functions for Cz in OH]): 



+ 



2(x-y)(l-z) 

/ 3 (1) (^) -9i\x,z) _ g2 ] (x,z) -g£\y,z) 
2(x-y)(l-z) 2{x-yY 
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P™ 1 (x,y)=P« L1 (x,y,l), 

P^ y) = ^)-fi%), 
x-y 

A.2.3 Box 

The loop functions for Cb u and Cz d in ([7]): 

B L u L \x) = (*), 

B LL2 (x v z) - 2 [ fWfx) _ /^(aQ-gPW) 

_ - ^(z, z) - (x Q y) 

ff ] ( x ) - 9i (x, z) g { 2 ] (x, z) - g$> (y, z) 



B* L1 (x,y,z) 



3(x — y)\ 1 — z x — y 



B*»{x,v) = ~^-A f?\ X ) - -^0(*) - /2( ° )(X) - 



3(1 -y) |/° v 7 l-j/ V ' ' 1-2/ 



= i/ 6 (1) (x), 
^ L2 (^y) = -^ 7T 1 - T i^ 1) (x) 



3(1-2/) l JS v ~ 7 1-2/ 
A. 2. 4 Magnetic and chromo- magnetic coefficients 

The loop functions for CV 7 and C$g in 

D L E L1 (x,y) = —(f 7 (x)-f 7 (y)), 
x — y 

Df\ x) = -/«(*) + ^/f) w + M£), 

n jiE,i/ x _ Mg) - M^vg) - (s y) fea(g) ~ M^*) ~ (x y) 
£ 3(x-y)(l-*) + 2(s-y)(l-s) 
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1 h^)- f * ){x) - f > 1){x > z) -(^y)l 



(x-y)(l-z) [ J1 v ' 1-z 



,RL3,„ .a _ 1 \vt„\ h 2 (x) -h 2 (x,y) 



D^(x,y) = — -th 5 (x) 



3(1 -y) I 1-y 

1 ft / x &3 (z) - h*{x,y) 
Hq[x) 



2(1 -y) I 1-Z/ 

I^foy) = —(/■(*) -/-(y)), 
x — y 

D L G L2 (x,y) = - 3 -fi 1 \x) + lfP(x), 

^RLl/„ .. _"\ _ ^0*0 ~ ^O^,*) - (x -B- y) 



D% L1 (x,y,z) 



2(x-y)(l-z) 



(x-y)(l-z) [ 1-z 



■\RL2> / ^ „.\ 1 \ut~\ h 2 (x) -h 2 (x,y) 



where fr(x), fs(x) are given by 



. , , 13 - 7a; 3 + 2x - 2x 2 , 

MX) = ; ttt : -— log£, 

1/71 7 24(x- l) 3 12(x-l) 4 & ' 
, , , 1 + 5x x(2 + x) , 

/s( " ) = -^^^ + 4(^Ir l0gX - 



The loop functions for CV 7 and Csg hi (E 



3(1-2/) I 1-2/ 

+ 20^)V 7 (x) 1^ 

^) = -|/l (0) (*) + /i 1) (^)-^/7 {1) W, 



2(1 -y) I l-y 

^) = -A (0) W + ^4 (1) W. 
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